Geometrical indications of adsorbed hydrogen atoms on graphite producing starhke 
and elhpsoidal features in scanning tunnehng microscopy images 
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Recent scanning tunneling spectroscopy (STM) experiments display images with star and ellip- 
soidal like features resulting from unique geometrical arrangements of a few adsorbed hydrogen 
atoms on graphite. Based on first-principles STM simulations, we propose a new model with three 
hydrogen atoms adsorbed on the graphene sheet in the shape of an equilateral triangle with a 
hexagon ring surrounded inside. The model reproduces the experimentally observed starlike STM 
patterns. Additionally, we confirm that an ortho- hydrogen pair is the configuration corresponding 
to the ellipsoidal images. These calculations reveal that when the hydrogen pairs are in the same 
orientation, they are energetically more stable. 



Hydrogen adsorption on graphite surface has been the 
subject of many theoretical and experimental studies, 
due to its fundamental importance in both science and 
technology [J, i, i, i, i, & 0, H ■ The STM technique is 
known as one of the most powerful experimental tools for 
analyzing surfaces. Through this technique, for example, 
Hornekaer et al. considered the atomic structures of var- 
ious types of H(D) clusters on graphite [1, 0|- In their 
latest experiments, they identified two types of clusters as 
the most abundant hydrogen species on graphite surface. 
The STM images of these clusters exhibit star and ellip- 
soidal like features, see FIG.[lI^a) 8]. The former includes 
six bright spots, three of which are relatively larger than 
the others Q]. Considering the STM images observed at 
various bias voltages, Hornekaer et al. concluded that the 
starlike feature originates from a particular arrangement 
of three or four hydrogen atoms adsorbed on graphite 
surface [SJ- For the ellipsoidal STM feature, using a 
set of STM simulations, Hornekaer et al. proposed that 
it should represent a pair of hydrogen atoms adsorbed 
on two adjacent carbon atoms, so-called ortho-hydrogen 
pair. Although, the computed STM image by Hornekaer 
et al. appears to be in agreement with experiment, the 
authors have not clarified how a pair of hydrogen atoms 
with H-H distance ^2.1 A can produce an elongated- 
ellipsoidal feature with a length of ~7.0 A 1(3. . In 
the present work we propose a model whose computed 
STM image perfectly matches with experimentally ob- 
served star like STM images. Additionally, we confirm 
that the ortho-hydrogen pair is the configuration for the 
ellipsoidal images and explain how these pairs produced 
the elongated features in STM images. We further in- 
vestigate the interaction of hydrogen pairs together. Our 
calculations show that when the pairs are oriented in the 
same direction, they are energetically more stable. 

The electronic structure calculations are carried out 
within the context of density functional theory (DFT) us- 
ing the spin-polarized Perdew— Burke— Ernzerhof (PBE) 
exchange-correlation functionals and the projected aug- 
mented wave (PAW) method, as implemented in the 
VASP code 0. The adsorbed hydrogen atoms on 
graphite surfaces are simulated using a graphene sheet 



containing 96 carbon and appropriate number of hydro- 
gen atoms (two, three, or four). The corresponding Bril- 
louin Zone is sampled by a 5 x 5 x 1 Monkhorst-Pack mesh. 
The structures are fully optimized until the magnitude of 
force on each ion becomes less than 0.04 eV/A. There is 
currently no simple way to apply the electric field in the 
VASP calculations. Accordingly, we use the SIESTA, 
which is a DFT code with localized basis sets for 
the simulation of STM images. We performed a new set 
of spin-polarized single-point energy calculations for the 
structures, previously optimized by VASP, in the same 
level of theory (DFT/PBE) but at definite electric fields. 
A full description of our methodology for STM calcula- 
tions can be found in Ref. pTllT^fTsj. 

To identify the geometrical positions of hydrogen 
atoms on graphite producing the starlike STM image, 
we first examine models with four H atoms (labeled as SI 
and S2 ) and three H atoms (labeled as S3 and S4), pro- 
posed by Hornekaer et al. 0], see FIG.[lJb). Figures[ljc- 
f) show the STM images computed for SI, S2, S3, and S4 
under an electric field of 0.5 V/A. From FIGs. [Ijc-f), it 
is evident that none of the above models can reproduce 
the experimentally observed STM pattern. Interestingly, 
the computed STM images of SI and S2 with four hy- 
drogen atoms have only three bright spots while the S3 
structure with three hydrogen atoms shows four bright 
spots in its STM pattern. This can be attributed to the 
fact that the computed current is substantially propor- 
tional to the electron tunneling and LDOS. For SI and 
S2, our calculations show almost no LDOS at left turn- 
ing points of the central hydrogen atom. As a result, 
the emitted current from the central hydrogen atom be- 
comes negligibly small. On the other hand, for S3, both 
the electron tunneling and LDOS turn out to be signifi- 
cantly higher for H atoms and their central carbon atom. 
Accordingly, these four atoms contribute substantially to 
the total emitted current, as can be seen in FIG. [He). 

We propose a new model S5, in which three hydrogen 
atoms are symmetrically placed on the graphene sheet in 
an equilateral triangle encompassing a complete hexagon 
ring of carbon atoms, see FIG. Wis)- The configuration 
S5 has a large magnetic moment of S/i^. 



(c) (d) 



• 


(0 










(h) 










•••• 












(n 








Ci) 


• • 
• • 

• • 


• 

• 

• 



FIG. 1; (Color online) (a) Experimental STM image of ad- 
sorbed hydrogen atoms on graphite with starlike and ellip- 
soidal features §]. The magnified starlike feature is shown 
in the inset, (b) the experimentally proposed configurations 
for the adsorbed hydrogen atoms producing the starlike STM 
images Q, and (c-f) the computed STM images for S1-S4 
at 0.5 V/A, respectively. Red (blue) color denotes the high- 
est (lowest) intensity of the emitted currents, (g) Geometrical 
structure of S5, and (h-j) its corresponding STM images com- 
puted at 0.7, 0.5, and 0.3 V/A, respectively. 



The spin-polarized calculation shows that total energy 
is approximately 0.1 eV less than that estimated by non 
spin-polarized calculation. Therefore the S5 structure is 
expected to be magnetic even at room temperature. Fig- 
ures [TJh-j) show the STM images of S5 computed at 0.7, 
0.5, and 0.3 V/A, respectively. The computed STM im- 
ages all show a starlike feature similar to that observed 
in experiment. There are six bright spots, three of which 
look relatively larger than the other three. Comparing 
the position of the spots and the geometrical configura- 



FIG. 2: (Color online) The geometrical configuration of S6 
and its possible transformation pathway to S5. 



tion of S5, the centers of the larger (smaller) spots turn 
out to be on the hydrogen atoms (the carbon sites indi- 
cated by black spots in FIG.[Hg)). In agreement with ex- 
perimental results reported by Hornekaer et ai, our sim- 
ulations demonstrate that when the strength of applied 
voltages decreases (here from 0.7 to 0.3 V/A), the inten- 
sity of the small spots reduces. To explain the low inten- 
sity of small spots at low voltages, we have considered the 
patterns of total LDOS and total tunneling probability of 
85 (not shown here). These results show that the LDOS 
patterns are almost the same under various applied volt- 
ages. At high-applied voltage 0.7 V/A, the tunneling 
probabilities of electrons from both hydrogen and carbon 
atoms are large. At lower voltages, 0.3 and 0.5 V/A, the 
tunneling probability of electrons from hydrogen atoms 
is relatively larger than that from carbon atoms. Such 
a difference can be attributed to the fact that the ad- 
sorbed hydrogen atoms create dipole moments on the 
graphene surface whereby the effective potential around 
the hydrogen atoms decreases and consequently, the tun- 
neling probability of electrons from them increases. This 
explains the reason for high electron emission from hy- 
drogen atoms in all ranges of applied voltages. It is worth 
mentioning that further calculations using a larger unit 
cell with a side length of '^35A reveals no change in the 
computed STM images and magnetic properties of S5. 

To compare the stability of S5 with other structures, 
we summarize in Table HI the respective adsorption en- 
ergy values of hydrogen atoms for optimized S3, S4, S5 
and S6 configurations (the latter is illustrated in FIG.[2]). 
The adsorption energy is defined as Ea—Egraphene+nH ~ 

Egraphene "^Efj where -£/(^rap/iene-t-n// ; E f^^pji^j^Q and E}{ 

are the total energies obtained for a graphene sheet with 
n adsorbed hydrogen atoms, a perfect graphene sheet 
and an isolated H atom, respectively. Overall comparison 
shows that the most stable structures, S3 and S5, are en- 
ergetically isomeric. However, since the STM image of S3 
is totally different from that observed in experiment, S5 
appears to be practically preferable to S3. This implies 
that there should be a post-adsorption mechanism af- 
fecting the structure of initially adsorbed hydrogen clus- 
ters on graphite, otherwise, S3 and S5 should statistically 
have had the same population, in practice. Thus, the ex- 
perimental conditions, within which the above STM im- 
ages have been observed, need to be carefully taken into 
account. In the experiment, S5 is observed in samples 
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TABLE I: Adsorption energy (Ea) values obtained for the 
triatomic configurations, S3-S6, and the dimer configurations, 
D1-D7. 

structures F,a{eV) structures Ea (eV) 



S3 


-2.63 


Dl 


-2.8 


S4 


-2.61 


D2 


-1.65 


S5 


-2.62 


D3 


-2.75 


S6 


-2.43 


D4 


-1.59 






D5 


-2.18 






D6 


-1.67 






D7 


-1.83 



which are annealed to 525K after a heavy deposition of 
H(D) atoms on the surface. Since the utiHzed H(D) beam 
is extremely hot (1600-2200 K), hydrogen clusters, in var- 
ious configurations including S1-S6, are created 0, 
However, by annealing the sample, many of these clus- 
ters are evaporated from the surface or change to other 
configurations. A statistical study by Hornekaer et al. 
reveals that, the relative abundance of two specific H(D) 
species with starlike and ellipsoidal STM features dom- 
inantly increases, as the sample is further annealed to 
570 K [8|. The process, described above, suggests that 
after annealing the sample, many H(D) species initially 
deposited on the surface transform into S5 configuration. 
For such a transition, basically, species with similar sym- 
metry but higher energy such as S6 are preferred. 

A possible transition pathway from S6 to S5 is de- 
picted in FIG.m In this picture. The diffusion reactions 
S6^A, A^B, and B^S5 are assumed to pass through 
the barriers AE^ = E:^g^ - Eg^, AE^ = E:^^^ - 
and Ai?3 = E^g2 ~ whose saddle points are at the 
transition states, TSl, TS2, and TS3, respectively. Note 
that, above all E* terms include the zero-point energy 
correction. To find the correct transition paths and the 
corresponding energy barriers, we use the nudged elas- 
tic band method [15| . As a criterion for determining 
the saddle points, the phonon eigenmodes of each TS are 
examined so that they have one and just one imaginary 
frequency. Accordingly, AE^, AE2 and AE^, are ob- 
tained to be 0.48 eV, 0.84 eV and 0.86 eV, respectively. 
The maximum influence of vibrational zero-point energy 
correction on barriers is less than 0.12 eV. For the sake 
of comparison, we have also calculated the energy bar- 
rier for desorbing each of hydrogen atoms, supposed to 
be displaced through S6^A, A— >B, and B— >S5 reactions. 
The respective desorption energies are 0.67, 1.0, and 1.2 
eV. Thus by annealing, the hydrogen atoms in S6 pre- 
fer to diffuse to S5 through the process described above 
rather than to be desorbed from the graphene surface. 

To further asses the possibility of above reactions at 
different temperatures, we have estimated the rate of re- 
actions at 300K and 500K within the transition-state the- 
ory [l6|. In this approach, the rate of reaction is defined 

- r=M^^^M^ -P {Wr) where and .f^ 
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FIG. 3; (Color online) (a) the geometrical structure of Dl- 
D7, (b-h) the respective STM images computed for D1-D7. 
(i) The projection of STM image of Dl onto its geometrical 
configuration. 



are the eigenmodes of initial state (IS) and transition 
state (TS), and as pointed out earlier AE* is the en- 
ergy difference between TS and IS after including zero 
point energy correction. On this basis, the respective 
values of r at 300K and 500 K are expected to be S6^A: 
11.2 X lOS-i and 3.1 x lO^s'^ A^B: 6.0 x IQ-^s'^ and 
4.3x10^5-1, andB^S5: 2.7xl0-^s-^, and 2.8 x lO'^s"!. 
These results clearly verify that the proposed transition 
process is practically very feasible. It is necessary to note 
that S5 may not be the most stable configuration for ad- 
sorption of three hydrogen atoms on graphene. Rather, 
it is a metastable structure made in a particular experi- 
mental conditions. 

As mentioned earlier, the experimentally observed 
STM images include patterns with elongated-ellipsoidal 
features [8]. Hornekaer et al. proposed that such fea- 
tures originate from a pair of H atoms, adsorbed on two 
adjacent carbon atoms of graphite @. Although their 
model, here labeled as Dl, appears to reproduce the 
same ellipsoidal STM features, their calculations fail to 
explain, how a pair of hydrogen atoms with H-H distance 
of ^2.1 A can produce a bright elhpsoidal feature as long 
as '^7.0A Q. To answer this question, we have computed 
the STM image for Dl and other possible H dimers. The 
structure of dimers and their corresponding STM images 
are illustrated in FIG.[3l Additionally, in Table[Tl we have 
summarized the values of adsorption energies. According 
to FIG. [21 each dimer has a unique STM pattern. Re- 
assuringly, Dl, is the only dimer structure whose STM 
image represents ellipsoidal feature, similar to that ob- 
served in experiment. Interestingly, the bright ellipsoidal 
pattern obtained in our STM calculation for Dl has al- 
most the same experimental length '^7.0A. For the sake 
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FIG. 4: (Color online) The relative positions of (a) two hy- 
drogen dimers and (b) three hydrogen dimers on a graphene 
sheet. The highlighted C-C bond denoted by numbers 1-9 
indicate the position of displaced dimer. 



of clarity, in FIG. [3i;i) we have projected the Dl's STM 
image onto its geometrical structure. Evidently, not only 
the hydrogen atoms but also their neighboring carbon 
atoms contribute significantly to the total STM current. 
In other words, the creation of Dl on graphene changes 
the electronic structure of its surrounding carbon atoms 
such that their corresponding LDOS values become sig- 
nificantly large. Consequently, they contribute substan- 
tially to the STM current. Dl is nonmagnetic system 
and according to Table HI it is energetically the most sta- 
ble configuration in comparison to other dimer models. 
Experimental results also show that Dl is the most abun- 
dantly formed configuration. 

Interestingly, it appears that many of the H dimers 
absorbed on graphite, seems to be oriented in the same 
direction. This can be clearly seen in FIG.[lJa), see the 
bright spots discriminated by the white circles. To elab- 
orate on this observation, we have carried out a set of 
calculations in which the adsorption energy is calculated 
for two and three H dimers when they have different ori- 
entations and distances in respect to each other. In our 
models as showir in FIGs.[l](a) and (b), we keep the po- 
sition of one or two of dimers fixed on the surface, while 
the other dimer is displaced so that for each calculation 
it is on one of C-C bonds, indicated by numbers 1-9. To 
minimize the interaction of dimers with their periodic im- 
ages in the neighboring unit cells, for calculations with 
two (three) H dimers we have considered a very large unit 



cell 34.08Ax 14.76 A (34.081x29. 52 A), containing 192 
(384) C atoms. The corresponding values of adsorption 
energies are summarized in Table [Til The results clearly 
indicate that dimers with same direction become energet- 
ically more stable, as they get closer to each other. On 
the other hand, the disorientation of H dimers results in 
an increase in the surface energy and, hence, in insta- 
bility of whole structure. Consequently, the adsorbed H 
dimers prefer to diffuse on graphene surface so that they 
can stay in the same directions. 

In conclusion based on STM image calculations, we 
have identified the geometrical configurations of the 
most abundant species of adsorbed hydrogen atoms on 
graphite after a heavy dosing of H(D) atoms. The struc- 

TABLE II: Adsorption energies Ea of hydrogen dimers on 
graphene sheets shown in FIGs. 21 a) and (b) 



Position of 


Ea of two dimers Ea 


of three dimers 


displaced dimer 


(eV) 


(eV) 


1 


-5.58 


-8.18 


2 


-5.49 


-8.21 


3 


-5.69 


-8.19 


4 


-5.54 


-8.18 


5 


-5.66 


-8.22 


6 


-5.39 


-8.16 


7 


-6.17 


-8.12 


8 


-5.64 


-8.27 


9 


-5.71 


-8.18 



tures were shown to have two and three hydrogen atoms 
with STM images having elongated-ellipsoidal and star- 
like features, respectively. The former (latter) turned out 
to be nonmagnetic (strongly magnetic). In the case of 
hydrogen pairs, they are energetically more stable when 
oriented in the same direction. 
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